Triplet contribution to the Josephson current in the nonequilibrium 
superconductor / ferromagnet / superconductor junction 



o 
o 

> 

o 



o 
o 

u 

o3 



i 

c 

o 
o 



> 



o 



X 



I. V. Bobkova and A. M. Bobkov 

Institute of Solid State Physics, Chernogolovka, Moscow reg., 142432 Russia 
(Dated: November 9, 2009) 

The Josephson current through a long s-wave superconductor/weak ferromagnet/s-wave supercon- 
ductor weak link is studied theoretically in the regime of nonequilibrium spin-dependent occupation 
of electron states in the ferromagnetic intelayer. While under the considered nonequilibrium condi- 
tion the standard supercurrent, carried by the singlet part of current-carrying density of states, is 
not modified, the additional supercurrent flowing via the triplet part of the current-carrying density 
of states appears. Depending on voltage, controlling the particular form of spin-dependent nonequi- 
librium in the interlayer, this additional current can enhance or reduce the usual current of the 
singlet component and also switch the junction between 0- and 7r-states. 

PACS numbers: 74.45.+C, 74.50.+r 



New states have been predicted and observed in 
Josephson weak links in the last years. One of them 
is the famous LOFF-state, which was predicted long 
ago [HQ- It was not observed in bulk materials, 
but can be induced in superconductor/weak ferromag- 
net /superconductor (SFS) Josephson junction, as it was 
predicted theoretically and observed experimentally 
[S SHU- The LOFF-state is qualitatively different from 
the ordinary zero-momentum superconducting state: as 
a response to the energy difference between the two spin 
directions Cooper pair acquires the total momentum 2Q 
or — 2Q inside the ferromagnet. Here Q oc h/vf, where h 
is an exchange energy and Vf is the Fermi velocity, repre- 
sents an increase (decrease )of the momentum of an elec- 
tron with the energetically favorable (unfavorable) spin. 
Combination of the two possibilities results in the spatial 
oscillations of the condensate wave function ty(x) in the 
ferromagnet along the direction normal to the SF inter- 
face. For the singlet Cooper pair ^ s (x) cx cos(2Q:c) @. 
The same picture is also valid in the diffusive limit. The 
only thing we need to add is an extra decay of the con- 
densate wave function due to scattering. In the regime 
h A, where A is a superconducting order parameter 
(OP) in the leads, the decay length is equal to the mag- 
netic coherence length £p = y/D/h, while the oscillation 
period is given by 2it^p. Here D is the diffusion constant 
in the ferromagnet, h = 1 throughout the paper. Because 
of the oscillations the condensate wave function contains 
nodes where the phase changes by ir. So, different signs 
of the OP can occur in the leads, what results in a nega- 
tive Josephson current response upon a small increase of 
the phase. This is the so-called 7r-state. 

However, manipulating by the phase difference x t» e ~ 
tween the OP's in the superconducting leads is not the 
only way to realize the 7r-state. The net supercurrent that 
flows between the two superconductors depends not only 
on the actual phase difference \ between the leads, but 
also on the occupation of the electron current-carrying 
states in the weak link. In the clean limit the current 



is carried by the Andreev bound states 13, 11, 12 1. In 
the diffusive limit the electron trajectories are not well 
defined, and Andreev bound states are no longer the 
appropriate concept to describe the supercurrent. The 
energy spectrum of the superconducting correlations is 
expressed in a so-called supercurrent- carrying density of 



states (SCDOS) [DjlHUJa]- This quantity represents the 



density of states weighted by a factor proportional to the 
current that each state carries in a certain direction. Pos- 
itive and negative parts of SCDOS give energy-dependent 
contributions to the supercurrent in the positive and neg- 
ative direction. The size and direction of the total su- 
percurrent depends therefore on the occupied fraction of 
these states, which is analogous to the occupation of the 
discrete Andreev bound states in a ballistic system. That 
is one can obtain negative Josephson current response 
to small phase differences and, hence, switch the system 
into 7r-state by creating an appropriate nonequilibrium 
quasiparticle distribution in the weak link region. This 



effect was predicted theoretically [13|, [1J] and observed 
experimentally [161 ] for a diffusive SNS junction. 

In the present paper we investigate the effect of the 
nonequilibrium occupation of the current-carrying states 
on the Josephson current in SFS junction in the param- 
eter range A <C h -C £/, where Ef is the Fermi energy 
of the ferromagnet. This regime is relevant to weak fer- 
romagnetic alloys, which were used for the experimental 
realization of magnetic 7r-j unctions. In this particular 
study we assume the Thouless energy stk — D/d 2 <C h, 
that is the interlayer length d ^> £p. It is found 
that a spin-independent nonequilibrium distribution of 
quasiparticles in the ferromagnet practically does not af- 
fect the Josephson current, that is the described above 
mechanism of the critical current reversal by the non- 
equilibrium redistribution of supercurrent-carrying states 
population is irrelevant in this case. However, the consid- 
ered system exhibits novel interesting physics upon cre- 
ating spin- dependent nonequilibrium quasiparticle distri- 
bution in the weak link. 
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It is well known that the presence of an exchange field 
leads to the formation of the triplet component of the 
condensate wave function in the interlayer. In the case 
of a homogeneous exchange field, which is considered in 
the present paper, only the component with zero spin 
projection on the field direction S z = is induced. Com- 
bining the two pairs with the total momenta 2Q and 
— 2Q into the triplet combination we see that the in- 
plane (S z — 0) triplet condensate wave function com- 
ponent ^t(x) oc sin(2Qa;), that is oscillates in space with 
the same period as the singlet one, but is shifted by ir/2 
with respect to it. Here we do not discuss the other triplet 
components with S z — ±1, which are typically induced 
in case of inhomogeneous magnetization (See [17j | and ref- 
erences therein). While their influence on the Josephson 
current under nonequilibrium conditions is undoubtedly 
of great interest, they do not exhibit damped oscillation 
behavior, which is crucial for us in the present study. 

The proximity-induced triplet component of the con- 
densate function is an odd- frequency quantity [I3| ■ Con- 
sequently, SCDOS generated by the in-plane triplet com- 
ponent is even function of quasiparticle energy. The su- 
percurrent can be expressed via the current-carrying den- 
sity as [m 



j oc / deNj{e)ip(e) 



(1) 



where e stands for the quasiparticle energy, Nj denotes 
SCDOS and ip(s) is the distribution function. Since 
<p(e) — tanhe/2T in equilibrium, it is easily seen that 
in-plane triplet component does not contribute to the 
supercurrent in this case. It is shown below that if the 
distribution function becomes nonequilibrium and spin- 
dependent, the supercurrent carried by the triplet com- 
ponent in the ferromagnet is non-zero. It should be noted 
that this current jtr does not compete with that one car- 
ried by the singlet component j s , but is additional to it. 
The magnitude of j tr can be of the same order or even 
larger than j s . Due to the fact that the singlet and triplet 
components have the same oscillation period but shifted 
in phase by tt/2, j tr can increase the usual supercurrent, 
carried by the singlet component, or weaken it, or even 
reverse the sign of the total supercurrent, thus switching 
between and 7r-state. Experimentally the most proba- 
ble way to realize the spin-dependent nonequilibrium in 
the interlayer is to apply a voltage to (or to pass the dis- 
sipative current through) a spin-active material. Then 
under appropriate conditions even quite small voltages 
should be enough to switch the system from to 7r-state 
and vice versa. 

For a quantitative analysis we use the formalism of 
quasiclassical Green-Keldysh functions in the diffusive 
limit 18, Il9j |. The fundamental quantity for diffusive 



transport is the momentum average of the quasiclassical 
Green's function g(x, e) — (g(pf,x,e)) Pf . It is a 8x8 ma- 
trix form in the product space of Keldysh, particle-hole 



and spin variables. Here x is the coordinate measured 
along the normal to the junction. 

The electric current should be calculated via Keldysh 
part of the quasiclassical Green's function. For the plane 
diffusive junction the corresponding expression reads as 
follows 



8Tr 2 eR F 



cfeTr4 



— (g{x,e)d x g{x,e)) K 



(2) 

where e is the electron charge and Rf stands for the re- 
sistance of the ferromagnetic region. (g(x , e)d x g(x, e)) K 
is a 4 x 4 Keldysh part of the corresponding combination 
of full Green's functions. r t and tTj are Pauli matrices in 
particle-hole and spin spaces respectively. 

It is convenient to express the Keldysh part of the full 
Green's function via the retarded and advanced compo- 
nents and the distribution function: g K = g R (p — (pg A . 
Here the argument (x, e) of all the functions is omitted for 
brevity. The distribution function is diagonal in particle- 
hole space: (p = <p(r + r 3 )/2 + g^gi { t o — T 3)/2- All the 
matrices denoted by " are 2x2 matrices in spin space 
throughout the paper. In terms of the distribution func- 
tion the current @ takes the form 



— oo 

-f R a 2 dJa 2 f A + (g R d x g R + f R dJ R )y - 
-<f{g A d x g A + f A dJ A ) 



(3) 



We assume that the direction of the exchange field h 
is spatially homogeneous and choose the quantization 
axis along the field. In this case the distribution func- 
tion and the normal part g R ' A of the Green's function 
are diagonal matrices in spin space. The anomalous 
Green's functions are represented as f R,A = f^' A i<J2 and 

f R ' A = — i<J2fd' , where f R ' A and f R,A are diagonal in 
spin space. That is the last two terms in Eq. ^ can be 

rewritten as (g R d x g R + f R dJ R )-g A d x g A - f A dJ A )£ cx 
Nj(e)tp. So these two terms are responsible for the su- 
percurrent carried by SCDOS Nj(s) occupied according 
to the distribution function ip(s). The first three terms 
originate from the direct diffusion of quasiparticles and 
are absent in equilibrium. 

The retarded and advanced Green's functions are ob- 
tained by solving the Usadel equations supplemented 
with the appropriate boundary conditions at SF inter- 
faces [2(|. We assume the junction to be long, that is 
rf>(p. This condition allows us to find the solution an- 
alytically even for an arbitrary SF interface transparency 
and low temperature, that is in the parameter region, 
where the equations cannot be linearized with respect 
to the anomalous Green's function. We start from the 
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completely incoherent junction (that is consider the left 
and right SF interfaces separately) and then calculate the 
corrections up to the first order of the small parameter 
exp [— d/£p] to the Green's functions. Within this accu- 
racy the anomalous Green's functions in the vicinity of 
left (1) and right (r) SF interfaces take the form 



Ida =™ 7T 



sinhe^ A e- Mx/2 +4S^ A (x)e 



iax/2 



J da 



~fda A (X -» -X) 



(4) 



Here a =t, J, is the electron spin projection on the quan- 
tization axis, k = +1(— 1) corresponds to the retarded 
(advanced) functions and a = +1(— 1) in the vicinity of 
the left (right) SF interface. The first term represents the 
anomalous Green's function at the ferromagnetic side of 
the isolated SF boundary with 



smhG*' A = 



4E% 4 (a:)(l 



(*)) 



(1-E# (*)) a 



(5) 



while the second one is the first order correction, orig- 
inated from the anomalous Green's function extended 
from the other SF interface. 

a a r> a — ( — \- (—)ax){\ + ina) / Ef 

X$,M= K ? ,A * 2 > (6) 

where K^ A is determined by the equation 

(1 + ina)K^ A {l - = |0 [sinhef^(l+ 

+K? Ai + Kf AA ) - coshGf' A 4K?> A (l + K?> a2 )] (7) 

R g stands for the resistance of each SF interface, which 
are supposed to be identical. sinhO^' and coshO^ ,A 
originated from the Green's functions in the supercon- 
ducting leads, which are assumed to be homogeneous and 
equilibrium. In this case 



R,A 



cosh0 



sinhO^- 4 



-me 



a/A 2 - (e + ni8) 2 
— niA 
a/A 2 - (i + KlS) 2 



(8) 



where 5 is a positive infinitesimal. 

The noncquilibrium distribution function in the inter- 
layer is supposed to be created by applying a voltage 
along the y direction between two additional electrodes 
N b and N t , which are attached to the central part of the 
interlayer. It is assumed that the conductances of N b F 
and N t F interfaces are spin-dependent and equal to g b(7 
and g t<7 , respectively. The voltage V ta (Ho-) between the 
superconducting leads and N t (N b ) electrode can also be 
spin-dependent. It can be realized, for example, by tak- 
ing the electrodes N b and Nt to be ferromagnetic and 
passing the dissipative current between them. Then the 



expression for the distribution function in the middle of 
the interlayer to the zero order of the parameter 
exp [—<1/£f] (disregarding the proximity effect) takes the 
form 



<P. 



(0) 



tanh ■ 



-eV tc 



2T 



gtaivF +d y g bo 



<?F(gta +9ba) 



+ tanh ■ 



-eV b . 



2T 



■9ba{°F +d y g ta ) 



-2d y g t(7 g b(7 



(9) 



where T is the temperature, o/ is the conductivity of F 
interlayer and d y is its width in the y direction. (p§' is 
connected to ip^ by the relation <p^](e) — — ipyt(— e). 

However, for the considered problem it is not enough 
to only calculate the distribution function without tak- 
ing into account the proximity effect in the interlayer. 
It should be obtain up to the first order of the param- 
eter exp [— d/^p]- For this purpose we solve the kinetic 
equation for the distribution function, which can be eas- 
ily derived from the Keldysh part of the Usadel equa- 
tion . The boundary conditions to the kinetic equation 
are also obtained from the Keldysh part of the general 
Kypriyanov-Lukichev boundary conditions [20| . 

Substituting the Green's functions ((4]) and the distri- 
bution functions into Eq. ^ we find the expression for 
Josephson current, which in the limit of low temperatures 
T <C |eVt,b | an d for |eX4 ,b| , | < A (the last condition al- 
lows us to avoid flowing of the dissipative current from 
the additional electrodes to the superconducting leads) 
takes the form 



J 



2eR a 



xIm(S^(-a-)sinhef) 



(10) 



Eq. ([10]) is the central result of our paper. While 
it is valid for arbitrary SF interface transparency, we 
concentrate on the discussion of the tunnel limit g = 
RF^F/Rgd -C 1, where Eq. (0 can be easily solved 
and the integral over energy can be calculated analyti- 
cally. For the analytical analysis we choose the most sim- 
ple form for the distribution function Eq. @ by setting 

g ta = and V b[ = 0. Then <^ 0) = tanh[(e - eV b ^)/2T], 

while ip'f — tanh[e/2T] and we obtain for the Josephson 
current the following result 



= i?f^Fsiny Ac _ d/fF 



AeR 2 g d 



i— . d 7T 

V27TC0S(- h —)- 

SF 4 



7i log 



A + eVfct 



eV b1 



. d 7T. 

Sn % + 4 } 



(11) 



The first term in Eq. (jlip represents the part of the su- 
percurrent carried by the singlet component of SCDOS. 
Under the conditions T <C |eVt,&| and |eVt,6^ J < A it 
is not affected by the fact that the distribution function 
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is nonequilibrium. The reason is that for a long junc- 
tion and h A the singlet part of SCDOS is concen- 
trated in the narrow energy intervals around e = ±A, 
where the distribution function does not differ from the 
equilibrium value. The second term j tr is caused by the 
triplet component and for qualitative understanding can 
be rewritten as 



Jtr oc 



*i 0) + ~(^ 0) 



■ (12) 



where the triplet component of SCDOS (Nj) tr ,a o 
d/dF _Al_ s i-a[d/iF +7r/4]. The first term in Eq. (TT2 



ere 



can be interpreted as the part of the supercurrent car- 
ried by the triplet component of SCDOS. Contrary to the 
singlet component of SCDOS, which is an odd function 
of energy, (Nj) tr ,^( e ) = ( N j)tr,a(-£) and, consequently, 
does not contribute to supercurrent in equilibrium. It is 
worth to note here that the spin-independent distribution 
function also does not give rise to the flow of the triplet 
component, because the total combination of the distri- 
bution functions in the square brakets is odd function of 
energy in this case. 

It is well known that the spin supercurrent can not 
flow through the singlet superconducting lead. On the 
other hand, it is seen that for our choice of the distribu- 
tion function the part of the supercurrent determined by 
the first term in Eq. (|12[) is carried entirely by spin-up 
electrons. However, the constant gradient of the distri- 
bution function oc (Nj) tr appears in the interlayer and 
cause the diffusive flow of electrons, which compensate 
the spin current. This flow originates from the first three 
terms in Eq.([3]) and is described by the second term in 
Eq. (H2D. 

It can be easily seen from Eq. that by manipulat- 
ing Vfe-f one can enhance or reduce the Josephson current, 
or to switch the junction between 0- and 7r-states. 
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FIG. 1. The critical Josephson current as a function of 
eV&j/A. The current is measured in units j c \e\R g / A. The 
different curves correspond to different lengths d of the junc- 
tion, which are measured in units of £f- 

To calculate the Josephson current for arbitrary trans- 
parency of SF interfaces one needs to solve Eq. ([7]) nu- 
merically and make use of Eq. (|10[) . The corresponding 
curves are plotted in Fig. [TJ Panel (a) shows the cur- 
rent for low enough dimensionless conductance of SF in- 
terface g — 0.1, while panel (b) represents the case of 



highly transparent interface ej = 3.0. Different curves 
correspond to different length el of the junction. It is 
seen that the discussed above tunnel limit qualitatively 
captures the essential physics. In dependence on Vb^ the 
current can be enhanced or reduced with respect to its 
value at Vf>t = 0- If t ne length of the equilibrium junc- 
tion is not far from the 0-tt transition, then small enough 
voltage can switch between the states. 

In conclusion, we have studied the effects of nonequi- 
librium spin-dependent electron distribution in the 
weakly ferromagnetic interlayer on the Josephson current 
through SFS junction. It is shown that the nonequilib- 
rium spin-dependent electron distribution gives rise to 
the supercurrent carried by the triplet component of SC- 
DOS. Depending on voltage, controlling the particular 
form of spin-dependent nonequilibrium in the interlayer, 
this additional current can enhance or reduce the usual 
current of the singlet component and also switch the junc- 
tion between 0- and 7r-states. The principal difference of 
the considered effect from already discussed in the lit- 
erature voltage driven SNS transistor 3, 3 15, EH is 
that in our case not the population of the singlet part 
of SCDOS is redistributed leading to diminishing of the 
current and reversing its sign, but the additional triplet 
component begins to contribute to supercurrent. 
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